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THE CACUULATIOH OF DRAG FOR AIRFOIL SECTIONS AHD 
BODIES OF EEYOLUnOH AT SUBCEITICAL SPEEDS 
By Max. A. Heaslet and. Gerald. E. Nitzberg 


SUMMARY 

A method, is developed for calc tile t ' jg the dreg, in a real 
compressible fluid and at subcritical Mach numbers, of airfoil 
sections at arbitrary lift coefficients and of bodies of revolu- 
tion at zero angle of attack. To apply the method it is necessary 
to know the velocity distribution for airfoils and the velocity and 
thickness distributions for bodies of revolution, together with the 
Mach number of the free— stream transition point from laminar to 
turbulent flow, and the Reynolds number based on chord or axial 
length. The method consists of tracing the growth of momentum 
thickness along the surface, for both the laminar and turbulent 
boundary layers, by means of relations which involve elementary 
integrals and can be evaluated by simple numerical means. An 
outline of the computational procedures required for drag calcula- 
tions is presented in the appendix to the report. 

The values of drag coefficient, computed hy the method of the 
present report for a number of cases, are compared with the values 
obtained for the same configurations by other methods and the dif- 
ferences between the various results are found to lie within the 
limits of accuracy of current experimental techniques. ‘ The use of 
the present method is recommended by its simplicity and generality. 


: INTRODUCTION 

Starting with the work of Prandtl (reference l), which was 
designed to determine the skin friction on a pointed flat plate in 
a uniform incompressible two-dimensional flow, the theory of drag 
calculations has been extended by several Investigators so that, 
under controlled “conditions and at speeds where air may be assumed 
an incompressible medium, very good agreement has been obtained 
with experiment for both airfoil sections and streamlined bodies 
of revolution. The calculation of drag is, however, limited to 
cases for which it is possible to estimate the location of the 
transition point, that ie, the point at which the laminar boundary 
layer over the forward portion of the body is terminated by the 
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onset of turbulent flow, and for whioh there is no extensive separa- 
tion of the turbulent boundary layer. 

In the present report the compressibility of the medium is 
considered and expressions for profile drag of airfoil sections 
and bodies of revolution at subcritical Mach numbers are given in 
forms which are particularly amenable to numerical calculation. 

The principal contribution, however, is contained in the treatment 
of the turbulent boundary layer in the two cases. As in previous 
work on this subject, the solutions consist essentially of inte- 
grals of the Kerman moment m equation for bodies in two-dimensional 
flow and for flow over three-dimensional bodies with axial symmetry. 
In reference 2, Squire and Young solve the problem for incompressible 
flow in two dimensions by means of a point -by— point method of 
integration requiring considerable labor, and In references 3, k, 

5, and 6 modifications of the Squire and Young method ore given in 
various forms which expedite the calculations. All these references 
give results which are in close agreement-. The method of Kalikhman 
in reference 6 is of particular interest for it is capable of 
generalization to the body of revolution and to the case of high- 
speed flow where density changes are of sufficient magnitude that 
they must be taken into account. This approach is adopted in the 
pre sent report . 

The various procedures which have been developed for predict- 
ing the growth of the turbulent boundary layer over an airfoil are 
all based on the same boundary— layer momentum aquation. In order 
to apply this equation it is necessary first to relate the skirt- 
friction coefficient to the boundary-layer momentum thickness. On 
the basis of experimental data for flat plates two such relationships 
have been evaluated: a power law (reference 7 ) and a logarithmic 

law (references 8 and 2). After comparison with the experimental 
data shown by Falkner In reference 'J , for Reynolds numbers between 
2 x 10 s and 5 X 10 7 , it appears that there is little significant 
difference in the numerical values of these two relations, when the 
scatter of the experimental data Is taken; Into consideration. The 
logarithmic law can be generalized easily to the case of compress- 
ible flow and Is used In the analysis of this report. 

The logarithmic relationship between the skin— friction .coeffi— . 
cient and the boundary— layer momentum thickness was combined by 
Squire and Young with the boundary— layer momentum equation to 
obtain the section drag of airfoils. The step— by— step Integration 
of the fundamental equation was first avoided in reference k where 
it was found that a considerable simplification can be achieved by 
dividing the velocity distribution over the airfoil Into segments 
In each of which the chordwiso velocity gradient Is relatively 
constant. Then, using an average value of the velocity gradient 
for each segment, it was found possible to construct a general graph 
from which the solution for any velocity distribution can be read. 

The authors of the present report were able to generalize the method 
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of reference if- to the case of compressible flov over airfoil sections 
but the" results have never been’ published. Xt was thought that the 
closed form in which the present results are given* together with 
the duality which it was possible to establish between the two— and 
three— dimensional cases, make "the latter approach .preferable. 

The power law relationship postulated by Falkner has been used 
in reference 3 by Holt to obtain a directly integrable relation for 
the turbulent boundary— layer growth. By means of a theoretical 
approach based on experimental results Tetervin (reference 5) 
related the skin-friction coefficient and boundary— layer momentum 
thickness in a more complex form which varied with the boundary- 
layer Reynolds number. Approximating this expression, over the 
range of integration, by a power law, Tetervin was able to express 
the growth of the turbulent layer in a manner somewhat analogous 
to that of Holt. The final forms resulting from this method of 
approach share with the present results for the turbulent layer the 
advantage of being in closed form. In reference 9 Tetervin has 
extended his method to include both two— and three-dimensional 
compressible flow. 

In the vicinity of the airfoil leading edge there is always 
a more or leBS extensive region of laminar boundary— layer flow. 

For airfoils at flight Reynolds numbers the laminar portion of the 
boundary layer contributes a minor portion of the total section drag; 
however, the amount is usually not negligible. In reference 10, 

Young and Winterbottom present a method for laminar boundary— layer 
calculations which includes compressibility effects. The derivation 
of their method is comparable to that of reference 11. There are, 
however, two signif leant differences: First, roference 10 Is based 

on Pohlhausen* s relationship for the velocity variation through the 
boundary layer, while reference 11 uses the Blasius velocity profile; 
and, Becond, reference 10 neglects the fact that for air Prandtl ’ a 
number la not equal to unity: The method of reference 11 Is used 

in the present roper t. 

Most of the theoretical and experimental work on bodies of 
revolution to date has been on airship shapes. With the present 
trend., however, toward large land-hased airplanes, particularly 
those with pressurized cabins. It is to be expected that fuselage 
shapes will approach bodies of revolution. The problem of studying 
the boundary— layer growth and the drag of bodies of revolution thus 
takes oh increased significance while at the same time It becomes 
necessary to generalize the procedure to include the effects of 
compressibility. The development of the' land. nary boundary layer 
over bodies of revolution in a compressible fluid Is given In 
reference 11 and the theory given there is applied directly In the 
present report. The momentum equation of the turbulent boundary 
layer Is given by Yeung in reference 12 for zero angles of incidence 
end a stop— by— step method of integration Is presented whereby the 
growth of the boundary layer may be determined for incompressible \ '■ 



k 


KACA EM No. A7B06 


flow. The boundary— layer equations, for "both turbulent and laminar 
flow, eire more complicated, for the body of revolution than for on 
airfoil section because of the fact that is is necessary to take 
into consideration the variation of the body radius along the axis. 

Drag calculations for bodies -of revolution have not been 
studied as extensively aB for airfoils and little previous work in 
the field of compressible flow has as yet been published. The 
present theory is similar to that developed for airfoil sections 
in that momentum loss in the boundary layer is expressed as a 
definite integral but differs in that it becomes necessary to modify 
the theory over the far aft portion of the body. In spite of this 
difficulty the method given does curtail sharply the amount of time 
required for the total calculation. . ■ _ • ■ 

A complete list of symbols, as used throughout tHlsreport, 
may be found in Appendix A, and tho computational- procedure for 
drag calculations is presented in Appendix B. , 


- "THEORY . ' . • 

Airfoil Sections 

Introductory remarks .— In figure 1 the two-dimensional flow 
about an airfoil section is indicated along with the boundary layer 
and wake associated with the flow. It is an established practice, 
in all theory connected with the calculation of drag, to divide the 
boundary layer and wake region into three different regimes of flow. 
Thus, if S represents tho stagnation point, tho boundary layer 
between S and the transition point at T.P, on either surface is 
laminar while between T.P. and the trailing edge at T.E, a turbulent 
boundary layer exists. In tho wake, the third region to be 
considered, the plana AA is drawn normal to tho center line of the 
wake at the point where static pressure in the wake has returned to 
its original free— stream value. 

It is easy to show, from momentum considerations, that if 
static pressure is assumed constant across tho wake, then tho drag 
D per unit length of tho airfoil is given by 

D = / w pu(U 0 - u)dy (1) 

whoro tho integration extends across tho wake in plane AA and 
u local velocity in wake 
p density in wako 
U 0 velocity of undisturbed stream 
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y distance measured normal .to center line of wake 

Momentum thickness of the boundary layer is, by definition. 



where 

U Telocity at edge of boundary layer 
u _ local velocity in boundary layer 
y distance measured normal to surface 
9 momentum thickness of boundary layer 
S boundary-layer thickness 
py density corresponding to velocity U 

and in a similar manner the momentum thickness of the wake may be 
defined. How let 

^-^ C 1 ^) 47 (3> 

where the integration is in plane AA and c 0 is density corresponding 
to free— stream velocity U 0 . Since drag coefficient c d is fixed 
by the relation 

D=c d |p 0 U 0 2 c (b) 


where c is. the chord length of the airfoil, it follows that 


c 


d 



(5) 


The analysis consists essentially in tracing the growth of 8 , 
the momentum thickness, along the top and bottom surface of the 
airfoil and in the wake to the plane AA. Since the nature of the 
flow in the boundary layer affects the rate of growth of the 
momentum thickness, it is necessary to treat the different regimes 
separately. The following development is therefore arranged to 
conform with this natural division. 

Laminar laye r.— In reference 11, expressions have been developed 
which may be applied immediately to determine the growth of the 
laminar boundary layer in two-dimensional flow of a compressible 
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fluid. In tliiB reference tha b oundary— laye r thickness di at the 
point Xi is defined as the distance from, the surface of the 
airfoil to a point in the "boundary layer where the ratio of the 
local, velocity to the velocity outeide the "boundary layer is 0.707. 
Neglecting terms involving the fourth power of M, where M is 
the Mach number of the free stream, the boundary— layer thickness 
dj. is given by the relation 


GKV 




(^) 


9 • IT 


5.3|^.35m{i-1.9i(^) 2 ]| 


4 (§)l Gg-rX#) -°- w (£T* T <?)" 


( 6 a) 


where , . ... 

E c Reynolds number based on chord length 
TJ a velocity outside boundary layer at point Xi 
U velocity outside boundary" layer at point x 
x distance along airfoil chord 


In the computation of section drag coefficients, inasmuch as 
the laminar portion of the boundary layer contributes a minor 
portion of "the total drag, it is practicable to simplify this 
equation. The modification will concern itself with the last term 
in equation ( 6 a) and is Justified by the fact that the last term 
contributes a small part to the tota.l value of. the boundary— layer 
thickness squared. Approximating the last term by the expression 
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Again neglecting terms involving the fourth povor of M, it is 
possible to show from results given in the reference that if 

Prandtl’B number Pr * is set equal to 0 . 733 * its value for 

free air* then 1 


(I ai)‘ - 0.082 ji + o.a «■ [1 - 1.3 j 


In the calculations that follow the basic variable will bo the 
nondimonsional product of momentum thickness and density. The 
value of this variable at the point x = xx, independent of the 
definition of boundary— layer thickness* is an immediate consequence 
of equations (6b) and (7)> Thus 


/ 0 PC\ 2 
^ c p o) x 


2i.l*aL ! 1 + 0.26 M 2 




dtx/c 

o \ u o/ 


Turbulent lay er.— The momentum equation for the turbulent 
boundary layer in compressible flow is given in reference 10 


JSL +f(H, 2 ) si * piffle „ -3-2 
ax 1 o' oa J ■ suit 


where the primes indicate differentiation with respect to x* 

H is a function of the boundary-layer velocity— profile shape* and 
t is the skin friction per unit area.. Under the assumption that* 
for compres sible fluids* 

* 0.2k5k o C * 3914 ^ ( 10 ) 

tV 

viiere 

f 2 = M! 

> T 


an& ix^ is the coefficient of viscosity at the -wall, it is 


1 
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possible to transform the equation to the f.-rr* "cod by Young and 
Winterbottom . ..T 


dx 


+ 2.555 (5 + 


1)21 = 22 lo.imi 

U 


-2 - 0 . 3914 ^ 

e 


( 11 ) 


The numerical mothods used in the integration of equation 
(ll) are somewhat protracted. To obviato this introduce now the 
transformation 


HPpy 

■ 


0 . 245 H e 


2 o » 39 ia £ 


( 12 ) 


With this change of variable, equation (ll) becomes 



+ K (K+l) 


[d (U/g 0 )/d(x /o)3 

U/U 0 


_ x- u UoPoc Pu 

u 


(13) 


where 

K = 2.555 (c.^914 + |) 

and the equation has been written in nondimensional form. 

It is necessary in equation (13) to relate the value of coef- 
ficient of viscosity at the wall to its value in the free stream 
M 0 . This follows ..directly from the empirical relation (re ference 13) 

*w VT W J 

whore T Q and T w are absolute temperature s in the stream and at 
the wall, together with 


■> *£- « 2 ) 

which is an immediate consequence of the assumption that energy is 
oonstant through the turbulent boundary layer. It follows that, 
approximately. 
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Setting 



and using the approximation for coefficient of viscosity, equation 
( 13 ) appears in final form 


+ K (H+l) I g n / d x l z = K U - P _— - -2 ( 15 ) 

41 U [l + 0.152 M 2 ] 

It is possible to put the solution of equation (15) in a form 
which is well adapted to calculations- if constant average values 
are used for E and K. Under this assumption, the integral of 
the differential equation is 


- if io * _Jk 


[ 1 + 0.152 M 2 J 


Ik 


_ _ K(H+x)+l _ 


p n 


dx 


} 


( 16 ) 


The variable in the - turbulent region, lies roughly between 

20 and 30 so that the total variation of JC Is small and K has 
an average value approximately equal to 1,21. The shape factor 
H varies, for a nonseparatod boundary layer, approximately from 
1.3 to about 1 . 7 , but from experience gained In other calculations 
It has been found that commutations. for low speeds are quite 
Insensitive to the value of E used and highly satisfactory 
results can be obtained for a constant value of E. In the present 
report K(H + l) shall be sot equal to 3- This assumes a value of 
E between 1.4 and 1.5 which is in conformity with low— speed 
measurements and, as shall be seen, will. give computed drags in 
close agreement with experiment and. other calculations. There are 
no available experimental measurements of velocity distributions 
through turbulent boundary layers at high speedB of sufficient ' • 
accuracy to permit the determination of the effect of compressibility 
on E. lacking such information, the assumption will be made that 
the same values for E can be used in the compressible case as In 
the incompressible case. Imposing the condition that at the transi- 
tion point z = Zj.p,* the arbitrary constant C Is determined and 

the solution becomes 



1.21 R r 


[ 1 + 0.152 


— P^- 


dx 


( 17 ) 


The density term In the integrand may be evaluated by assuming 
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tile flow outside the boundary layer is isentropic ; thus 


P 


1 + 2 =i M 2 (l-lj 2 ) 


1 

7-1 


(18) 


The value of z at the transition pointy which is required in 
equation ( 17 )* must he found from the value of determined 

by equation ( 8 ). Since* in non&imensional variables* 


where 


5 


2.555 


inf k -°75 R c g P g 'N 
\ [1+0.152 M 2 ]/ 


7.n deno tes natural logarithms, then 


z = 1.604 w ln®w 


(19) 


if 


w = 


ihon_ 

[1+0.152 M 2 ] 


E c 9 0 Tl 


( 20 ) 


Substituting from equation ( 8 ) into equation (20) thus gives v at 
the transition point, in figure 2 , which is a plot of equation ( 19 ) 
the required value of zy.p. can he found. If the velocity dlstri— 

hution over the airfoil and the free— stream Mach number are known 
it is now possible to substitute directly into equation ( 17 ) and 
determine the growth of the boundary layer up to' the trailing edge. 

Wake « — YounR and Winterbottom in reference 10 have jiiscussed 
the momentum equation in the wake and have concluded on the basis 
of what experimental data are available., that^ 

P2 02 = Pt.e. e T.E. U T.E. 3 * 2 ( 2i -* 


where subscript 2 applies at plane AA in the wake and T.E. indi- 
cates values at the. trailing edge of the airfoil._ 

From zipjj]. the value of w^g, follows and the airfoil 
section profile drag coefficient is given by 


c d = [1+0.1.52 M 2 ] r 


2wi 


T - E +- U. 


4.075 E, 


T.E. 


2.2 


( 22 ) 
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Bodies of Revolution 

Introdu ctory^ remarks.— - The de velopma nt of the theory relating 
to the calculation of drag for a Body of revolution is directly 
comparable to the theory presented above for airfoil sections except 
that the angle of attack qf the body will be restricted to zero. 

In figure 3 the body is shown; point S . representing the stagnation 
point, T.P. indicating the transition point in the piano of the 
paper, T.E. denoting t ho tail end of the body, and AA marking the 
position of the plane where static pressure has returned to its 
value in the ambient stream. 

Drag coefficient of the body is by definition 


o„ £ 

K D o *< T > e/3 

where 


D 

drag of body 



velocity in the free 

■ stream 

Po 

density in the free 

stream 

V 

volume of body 


and. 

from consideration 

of momentum. 


(23) 


D = 2 it/ pu(U 0 - u)y dy 

v 

whore the integration extends across the wake in plane AA and 
u local velocity in the wake 

p density in the wake 

y distance measured normal to center line of woke 

For bodies of revolution the momentum area 0 which is defined by 
the equation 

$ = 8 jc -f3£ ( 1 - (r + y ccs a)dy ( 25 ) 

J o 13 S 

where 

5 thickness of tho boundary layer 
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r radius of cross section of body 


a angle "between tangent to generator and axis of body 

occupies a role analogous to that of momentum thickness in two- 
dimensional airfoil theory. Momentum area in the wake, at the 
plane AA, is - 


whence 





(26) 


(27) 


The theory which follows will also have occasion to use the 
variable $ which is related to the momentum area by the 
expression 

Q = <t/2i tr (28) 


or 


* ^ C ~ § X 1 008 a ) 47 


(29) 


Laminar layer .— From, the theory developed in reference 11, the 
laminar boundary— layer thickness &i, defined as in the two- 
dimensional case, at an arbitrary point Xi is given, neglecting 
terms involving the fourth power of M, by the expression 




1 1-0.15 


M 2 



i 

/ (!) 


e d y 


:i7 d^S_o'.U M *J 


fEi 

l 


where 




Ui velocity outside boundary layer at point xi 
Rj Reynolds number based on length of body 


radius of cross section of body 


(30a) 


r 
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x distance measured along axis 

di "boundary— layer thickness as defined for equation. (6a) 

l axial length of "body 


Since, in general, the lanunar portion of the "boundary layer 
contributes a "small part of the total drag and- since the last term 
in equation (30a) represents a small part of the value given by 
the relation, it is practicable to derive a simplification for 
(di/i) s analogous to equation (6b). Thus, 



- 2-3 (. 

i\9* '■r V: 


E, 


A-Y* i I-O.35 ^ 

hi/ j 


[“•» ©-)■] 




(30b) 


It can alBO be shown that if Prandtl 1 s number is set equal 
to. 0.733 and if terms in M of fourth degree and higher are 
ignored, then for the’ tody of revol\rfcion the "following approximate 
relation is obtained. - r .-. • 


£ ^U\ z m 

d Po / i 


0.082 


1 


+- 0.6l M a 




( 31 ) 


Grouping momentum thickness and density together, sind in 
nondimensional form, the following equation holds 



Turbulent layer The momentum equation of the boundary 

layer, for a body of revolution in compressible flow is given 
in reference 12 in the form 
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T 2jtr 


where 

A* displacement area in boundary layer (see definition in list of 
symbols) 

t Bkin friction per unit of area 

Setting (A */<t) = H reduces the equation for the boundary layer of 
the bodies of revolution to a form similar to- that for airfoils. 
Thus 


dO 

dx 



P*U ~ 
PtI . 


$ = 


Ptfj 5 


5nr 


( 33 ) 


the primes indicating differentiation with respect to x. 


It should be noted that the definition of H differs in the 
three-dimenBional case from that in two dimensions. However, in 
the case where the thickness of the boundary layer is small in 
comparison with the local radius of the body of revolution the two 
expressions for H are approximately equal. 


The relation between t, p, U, and.„9 used in the analysis 
for airfoil sections was based on theory that held for a flat 
plate;.. that is, the pressure gradients we rs~' ignore d in that 
particular, phase of the study. Since on bodies of revolution the 
pressure gradients are small/ the same relation may be assumed to 
hold, thus 


Up£ 

lw 




O .3014t, 


m 


where 



From equation (3*0 and the definition of 0, it follows that 


Mw 


0.3914^ 

= 0.2k5k 2nr e 


( 35 ) 
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and substitution, of this relation into equation (33) gives 


+ 2.555 [ (H+l) ?- + El] 
d2 . J r J 


^Pg 

Pv 


=-2 — 0 . 3914 ^ 
10. hilt, e 


(36) 


Now let 


z 


110 PU f 2 
Pv & 


0.2lj-5^ 2 e° * 3914 ^ 


(37) 


Direct substitution into equation (36) yields 


-£U— + kF(H+ 1) + ft ( r / . l . )/d(r/l . ) 

d J - r/l 

U ^oPo 1 PU 


= K 


Po 


(38) 


vhero 


2.555 (0.391k + |^) = K 


V. 

and the equation is written in nondimens i onal form. 
Setting 


u _if pcj „ 7T x = v 
% ~ u * p ' T 


and using the relation 

= m 0 [1 + 0.152 

the final form of equation ( 38 ) "becomes 


M + K [(H+l) 

dx U r 


z • K 


UpR l 


[1+0.152 M 2 ] 


(39) 


Putting again K = 1.21, K(H+l) = 3 and fixing the arbitrary- 
constant of integration by assuming that z = at x = 

it is possible to express the solution of equation (39) in the 
form 
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z 




1.81 E; 

U 3 [l+0.152 m 2 ]? 1 - 21 


/ 


p U 4 r 1 * 21 dx 


^T.P. 


m 


The density term in the integrand is evaluated "by means of the 
expression 


P = 


1 + 7-1 M?(l 



(hi) 


The determination of z,p p from the known expression for 

(ffp) given hy equation (32). proceeds as follows: IVom equation 

(3M . 


i 


2.555 in / MN 

V [ 1 + 0.152 M 2 ] J 


whence 


z = 1.604 w In 


if 


= — _ e p u 
[1+0.152 M 2 ] 


(42) 

(43) 


Thus, from equations (t-3) and (32) the value of w at the transi- 
tion point can "be found and Zqi.p. is ohtainahle from figure 2. 

With this information, together with the velocity distribution 
over the body and the free—streaii Mach number equation (40) can be 
used to trace the growth of the boundary layer aft of the transi- 
tion point. ■ n. . ; ........... . 

One difference arises in the computations for bodies of revolu- 
tion which distinguishes the theory from that for airfoil sections. 
This is duo to the fact that r vanishes at the tail of the body 
and as a consequence an infinite singularity appears in equation 
(ho). Because. of this singularity it is not possible to carry the 
integration to the tail for momentum thickness’ will become infinitely 
large and the expression for drag coefficient becomes an indotermi— 
nate form. To circumvent this difficulty it is necessary to use 
equation (ho) up to some arbitrary point, say the 8o~percent point 
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of the axis, and then to modify the method of calculation. For 
this purpose it is convenient to compute the growth of momentum 
area rather than momentum thinknese over the latter portion of the 
turbulent run. 


Two means will he given whereby the momentum area can be 
calculated-- The first, which is merely an extension of reference 10 
to the. case of compressible flow, involves a point-by— point integra- 
tion of the basic differential equation. This equation is express-, 
ible nondimensi onally in the form 


dO 

dx 


dSr &P 

(H+2) 

IT p 


( 0.^914 ) a 2 it r 


In®- 


U p $ Ej 


m 


Q.'gkyh X 2rtr X [1+0.152 M 2 I 


I 


t 


The derivation of this expression may be obtained by combining 

a _ P£f ■ 

a T 


and 


Upn<& , , 0 . 39 J.4^ 

• —-P- = 0.2k5k 2rt r e 
Hv 


to g9-| 


( 0 . 391 k) ■ 


PU ? 2 - •- in 2 _ 


U p ® B z 


(k5) 


0.2k5k X 2 jt r X [1+0.152 M S I 


which, together with equation ( 33 )> will give the required relation. 

From the value of 6 and r at the 80-percent point on, the 
axis, the value of can be found at this point and the - growth of 
can then be calculated over the remaining portion of the body. 

In particular, if . is the value of $ at an arbitrary poinb 
~Z on the axis, • . 

r*n + (|| At*) , <«> 

where ^ n+ 2 is the value of ^ at the point x + A (2) . The 

calculation consists ■ of repeated applications of this relation. 

As the interval A(t) in equation (k 6 ) gets. smaller the result of 
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the calculation approaches the exact solution of equation (44). 

Since the ah ova calculation ie to be applied over an interval 
which is small in comparison with the axial length of the body, the 
labor of such' a calculation is much less than would be required if 
such methods were applied to the total turbulent run. It is possible 
however, to shorten this calculation further, by assuming that 
cf2rt r is a linear function of x where Cf = 2T/p 0 U 0 2 Is the 
local skin- friction coefficient. The validity of such an assumption 
will be examined later in the discussion. 


Return now to equation (33)* From the definition of skin- 
friction coefficient the right-hand side of this equation is 
expressible as 


nr 


2T Pp 
PoUo 2 pU 



(47) 


If Cf2« r is linear over the aft .portion of the body, falling 
from Its value at the 80-percent point to zero at the tail, then 
in nondimens ional terms equation (33) can be written as 


m + 

dx 


(H+2) (&/*=) + -(dp^dx) 


* =2 (°f 2rt? ) 0 . 8 


12*41 

pTT 


(48) 


This equation can be integrated, and as a result 

H 


_ _-(H+a) 

4 = U ■ 


P C +f | (c f 2n r) Q>8 (l-x)U dx 


(49) 


■When x = 0.8', 4 = ^ 0.8 

so that • ' • 

C = 


b l 1 

o.s 0.8 


H+2 _ 
P, 


0.8 


Using this value of C, together with equation (49), 

4 = 8 (%a) E+s ( £2^®)+ (i-x)Tj H dx 

o.eV TJ / V p J 2 ■* uH+2 J 


0.8 


(50) 


and, after substituting from equation (47), the value of 4 at the 
tail ie given "by the equation 
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TP. 35. 



3.5 


( 


Po.e \ 
PT.E.' 



■ 5 


2it r 

6.528 

l 

; 

^.075 Sj 6 p XT 

k. [1+0.152 1# ] 


) 


O *8 


5^,^_7 t -°(«)( = jL-y-- g .' ( 51 ) 

PT.E.^T.B. 2 J 0>8 ' ,U T.E.' / 


Wake . - In reference 12 Young has integrated the momentum 
equation in the vaike applying the same methods used in references 
2 and 10, for the body of revolution. If subscript 2 indicates 
values of the variables In the plane AA and TJS. denotes values 
at the tail of the body, then It is ehown that 


c 2 $2 = Pt.e. .^D.E. u t.e. 


3,2 


(52) 


Since, from' equation (2?) 


Cd = 


2®, 


2/3 


It follows that 


Cjj «. 


2 Pr.E. j T.E. 

Y a/a /i*; 



3.2 


(53) 


DISCUSSION 

The foregoing theory provides a convenient procedure for 
studying the growth of laminar and turbulent boundary layers and; 
for calculating the drag coefficients of -airfoils arid bodies of 
revolution. Just as in the case of wind-tunnel testing, where it 
is essential that the model tested be an accurate representation of 
the original configuration, it is Important that in the application 
of this. theory the operator should be able to determine correctly 
the required ‘aerodynamic properties of the configuration under 
consideration. * This implies that the pressure distribution over 
the body and the extent of the laminar* and turbulent layers be 
specif led or be determinable. 
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Telocity distributions are an immediate consequence of pressure 
distributions, for both low- and high-speed flow, "so that if experi- 
mental data are . aval lable_ the calculations may proceed directly. 

The theoretical calculation of the yeicfcity distribution corre- 
sponding to -a given shape is, on the other hand, a rather lengthy 
process although such methods have been treated adequately in the 
literature. For an arbitrary a'irfoil section at any desired lift 
coefficient the velocity distribution for incompressible flow can 
be found by the methods of Theodoreen (reference it) , Allen 
(reference 15 ) or Goldstein (reference 16) . For NACA conventional 
and low-drag airfoils corresponding distributions may be found quite 
easily from the tabular data given in reference 17* At subcritical 
Mach numbers the velocity distributions are calculable from low- 
speed data by moans of the well-known Glauert— Prandtl or Karman- 
Tsien transformations. For the body of revolution, methods have 
been given by Young and Owen (reference l6) and Kaplan (reference 19 ) • 

The theory for the determination of velocity distributions is 
of course based on the assumption of potential flow but, for 
airfoils, if the lift coefficient rather than angle of attack is 
specified the calculations are sufficiently accurate for most appli- 
cations. In reference 20 it is shown that the effect of the presence 
of a boundary layer is primarily to change the apparent angle of 
attack of the airfoil and to increase the local volocitios in the 
vicinity of the trailing ed®3. A procedure is introduced in this 
reference for estimating the magnitude of the change in the trailing- 
e&ge velocity brought about by the presence of the boundary layer. 

This procedure first estimates boundary— layer thickness from the 
potential theory velocity distribution and can be used in conjunc- 
tion with the theory of the present report. In the calculation of 
drag coefficient, however, the nature of the equation is such that 
the total drag coefficient computed is merely affected to a very 
small degree by moderate changes In the trailing-edge velocity, and 
as a consequence such a refinement is not used when only drag coeffi- 
cient is to be found. 

The determination of the location of the transition point from, 
laminar to turbulent flow in the boundary layer presents a diffi- 
cult problem. At small Reynolds numbers and for smooth surfaces 
transition occurs in a region of decreasing local volocitios where 
there is usually a region of separated laminar flow between the 
laminar and turbulent portions of the boundary layer. However, at 
Reynolds .numbers greater than several million the length of this 
region of separated laminar flow is of negligible extent so that It 
Is possible to consider transition as occurring at a point. Experi- 
mental flight tests of smooth airfoils with maximum velocity in the 
vicinity of the midchord Indicate that transition occurs when the 
local laminar boundary— layor Reynolds number R^ = Udp/p attains 
a value of about $000. The velocity TJ is the local velocity • 
outside the boundary layer and It should be noted in particular 
that the characteristic length d used in evaluating this Reynolds 
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number is the value of d obtained by the equation for the laminar 
boundary layer in the present report. This transition criterion 
applies strictly only to the determination of whether the transi- 
tion point is ahead of the maximum velocity point. For configura- 
tions in which maximum velocity is as far back as the midchord 
position and for which the transition Reynolds number occurs aft 
of this point, it is probable that the transition point will be 
close behind the maximum yelocity points For smooth airfoils at 
an g l es of attack such that there is a sharp velocity peak in the 
vicinity of the leading edge, transition, occurs behind the maximum 
velocity point. . Theory and experiment indicate that in such cases 
transition at large Reynolds numbers occurs only after the velocity 
has decreased between 5 and 10 percent of the maximum velocity, the 
percentage of the velocity recovery before transition occurs being 
greater the more slowly the velocity decreases in the chgrdvi se 
direction. 

The preceding criteria for the location of the transition 
point indicate the most rearward position that can -be expected, 
that -is, the probable position for smooth airfoils in low— turbulence 
flow. When the surface under consideration is rough or contains 
such transit.! on-promo t 1 ng agents as protuberances, waves, air 
leakage-, or dust particles ahead of the transition point, as 
predicted for ideal. conditions, it is to be anticipated that transi- 
tion will move forward in the direction of such, disturbances. In 
the immediate vicinity of the stagnation point, ;however, there is a 
very rapid acceleration of the air so that any local, disturbance 
which is not sufficiently severe to -change the local velocity’ 
distribution will be unable to cause transition to. occur in this 
region of very ' favorable velocity gradient. 

There are not sufficient experimental data on the lobation of 
the transition point' on bodies .of revolution in low— turbulence 
flows. Lacking such information, the most reasonable basis for 
estimating the transition point on bodies of revolution seems to 
be to use the same' criteria' 'as previously presented. for airfoils. 

Drag of 'Airfoil Sections 

In reference 2 the section dra^ coefficients, computed by the 
method of that report, are presented for an extensive range of 
airfoil thickness— chord ratios, Reynolds numbers, and transition 
point locations. A representative group of these cases has been 
recomputed by the methods of reference 4, as well as the present 
report, and the results of these calculations are given in table I. 
It is seen in the table that the maximum difference existing between 
the three methods for obtaining -the drag coefficient is 0.0002 for 
a single surface and is 0.0004 for the total drag coefficient for 
both airfoil surfaces. There is no apparent consistency in the 
nature of the deviations, however, and in only one case does the 
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difference in total drag reach the value given above. The results 
of reference 4 should agree numerically with those of reference 2, 
since they are based on the same fundamental assumptions, but some 
difference might he expected when comparisons are made vith the 
computations based on the present report because of the small 
change in the shape factor H and the averaging method used to 
fix the factor K. The agreement between the results should 
therefore be considered highly satisfactory and a confirmation of 
the compatibility of the assumptions. 

The limits of accuracy of. current methods for measuring 
airfoil section drag coefficients is of the same order of magnitude 
as the differences existing between the various theoretical results 
so that it is not possible to say which of the calculations most 
accurately predicts experimental values. The wake— survey method is 
now used commonly in the determination of experimental drag coeffi- 
cient and it can certainly not be said to determine drag within the 
limits needed to establish the relative accuracy of the preceding 
computations even though, for a given test configuration, it may 
be possible to repeat measurements to a higher order of accuracy. 

Any experimental, check i? also complicated by the problem of 
locating the point of transition from the laminar to turbulent flow 
in the boundary layer. The previously mentioned methods for 
determining the transition point can easily err by a few percent of 
the chord length on each surface, and this can bring about an error 
in the calculated section drag coefficient of the order of magni- 
tude of 0.0004. 

Very few experiments have been conducted in which the location 
of transition from laminar to turbulent boundary— layer flow on both 
surfaces and the corresponding section drag coefficient were measured 
accurately. In reference 2 the section drag coefficient 
measured in flight is given as O.OOSO for a 25— percent-thick 
section at a lift coefficient of 0.25 and a Reynolds number of 
8.2 x 10 s , the transition points having been measured and found to 
be 36 and ?,6 percent of the chord length, from the leading edge, 
on the upper and lower surfaces, respectively. For this configurar- 
tion the following results are - given by the _ various listed methods 
for calculating airfoil section drag coefficients: 

Method 'Section drag coefficient 


Experimental, Flight (reference 2) 

0.0080 

Squire and Young (reference 2) 

. 0079 

Nitzberg (reference 4) 

.0077 

Eolt (reference 3) 

.0076 

Tetervin (reference 5) 

.0077 

Present report 

.0080 



HACA EM No. A?B06 


23 


These value b are all within the probable limits of accuracy of 
experiment. . 

The drag coefficient for an NftCA 0012 airfoil at' zero angle 
. of attack and the corresponding transition-point locations . have 
"been measured at a. number of Reynolds numbers. These measurements 
are tabulated in reference 5 along with the theoretical drag 
coefficients obtained by several theoretical methods.' The results, 
together with the numerical values calculated as indicated in the 
present report, are px-esantod in table II.. It is believed that the 
differences between the numerical values of the various methods 
leave little choice as to the relative accuracy- Any decision to 
use one method in preference to the others must rest, for the 
present, on convenience of application. The procedure of the 
present report requires no approximations to the velocity distri- 
bution over the airfoil and it is readily applied to calculating 
both the growth of the laminar and turbulent boundary layers as 
well as computing airfoil section drag coefficient. 

The calculation of airfoil section drag coefficient for 
airfoils at speeds requiring the inclusion of compressibility 
effects was first given by Young and Vfi nterbo t tom (reference 10) . 

The single numerical example considered in this reference is an 
18 . 5 -percent— thick symmetrical Jbukowski section at zero angle of 
attack. The assumption was made that the transition point is 9-^ 
percent of. the- chord from the leading edge and that the Reynolds 
number is equal to 10 7 . Velocity distributions were used for 
potential flow and at a Mach number of 0.685- For these two cases 
the calculated drag coefficients were O.OO 89 and 0.0091, respectively, 
' while the present report gave, for the seme data, the values O.OO 89 
and 0.0093- Thus, both methods indicate that at subcritical Mach 
numbers the introduction of compressibility effects into the compu- 
tations brings about only a slight increase in the airfoil section 
drag coefficient. 


Drag of Bodies of • Revolution 

The problem of bound _,i‘y— layer growth and the determination- of 
drag coefficient for bodies of revolution is more complex than for 
airfoil sections, since it is necessary to take into consideration 
the dimensions of the body as well as the velocity distribution. 

The method derived in the present report furnishes a procedure 
which parallels closely the analysis derived for an airfoil section 
and, with little increase in intricacy, embraces compressibility 
effects. In order to compare results obtained by the present method 
with results given in reference If?, drag coefficients of the Akron 
airship shape were computed at a variety of Reynolds numbers and 
transition points; These results are presented in table HI along 
with the corresponding drag coefficients obtained by Young for the 
same configurations. (For convenience of comparison, -Young's 
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convention of "baaing the drag coefficient on the body surface area 
was adopted. This convention differs from those used by the'NACA 
which bases the drag coefficient oh either the volume to the two- 
thirds power or on the projected frontal area.) It is noted that 
for transition points far forward the difference between the results 
for identical configurations, as obtained by the present method and 
that of Young, increases with the Eeynolds number. This implies 
that the difference between the two methods arises in the calcula- 
tion for the turbulent layer and is probably brought about by the 
use of the averaging mat hod associated with the local skin-friction 
coefficient. 

Very few experiments have been conducted on bodies of revolu- 
tion at zero angle of attack to determine the location of transi- 
tion and the corresponding drag coefficient; In reference 21 the 
drag— coefficient measurements for the Akcpn airship shape at three 
Reynolds numbers of the order' of 10 , 000,000 are given. Boundary- 
layer surveys indicated that in each of these three cases the 
transition from laminar to turbulent boundary— layer flow occurred 
at 7 percent of the body length from the leading edge. In 
reference 22 drag coefficients are given at the same three Reynolds 
numbers for a metal body of virtually the same shape. These latter 
measurements are considerably larger, a fact which is difficult to 
explain because the pressure distribution over the forward 5 percent 
of the body is so very favorable that at these moderate test Reynolds 
numbers it is improbable that transition could move significantly 
ahead of the 7— percent station. For the series of measurements 
reference 22 indicates that the wind— tunnel— interference effect was 
of minor importance. The wooden model was of polygonal cross 
section; whereas the metal model was a true body' of revolution but 
this would seem to bo unimportant. Assuming' that transition 
occurred at 7 percent of the body length from the leading edge and 
using the experimental pressure distribution of reference 21 , the 
drag coefficient of the Akron shape was calculated at' the three 
test Reynolds numbers. 

In the following table the calculated values are compared with 
the two sets of experimental values. 


Reynolds 

number 

12.3 X. 1 Q 6 

15.0 X 10 s 

17.3 X 10 s 

Wooden 

model 

•0.0198 

0.0193 

0.0190 

Metal 

model 

.0228 

.0223 , 

.0219 

Theory 

.0222 

. . 0217 

! . . . - — j 

.0212 
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It is seen that the calculated values lie "between the two sets of 
experimental data. 

Raw data from wind-tunnel experiments on "bodies of revolution 
generally indicate a change in. the absolute value of the pressure 
coefficient as the Mach number increases but the magnitude -of this 
increase is such that it may very well be attributable to wind- 
tunnel— wall effects. There is a difference of opinion between 
various authors as to the effects of compressibility on the velocity 
distribution over bodies of revolution for analyses have been 
presented both affirming and denying that pressure coefficients 
rise with increasing Mach numbers. These differences have been 
resolved by J. G-. Harriot. In reference 23 he shows by means of 
linear perturbation theory that, for very slender streamline bodies 
of revolution in a uniform stream of compressible fluid, the' 
pressure coefficient at the surface of the body is almost independ- 
ent of Mach number. Hie equations of the present report, together 
with the result that the velocity distribution is independent of 
Mach number, were used to determine the effect of compressibility 
on the calculated drag coefficient for the Akron shaped For transi- 
tion point at 25 .7 percent of the chord from the leading edge, 
Reynolds number equal to 10 a , and free— stream Mach number equal to 
0.7, this configuration had for calculated drag coefficient the 
value 0.00184 as compared with O.OOI 98 for the incompressible case. 

In the presentation of the theory it was noted that the method 
developed for computing the growth of the boundary layer over bodies 
of revolution breaks down in the vicinity of the tail end. When 
the drag of the Akron shape was being computed it was observed that, 
over a range of Reynolds numbers from 10 s to 10 s and transition 
points from the leading edge to about the mid chord position, in the 
step— by— step integration over the last 20 percent of the chord the 
quantity c f [2ir(r/z}] varied almost linearly with x from Its 
value at the 80 --percent— chord station to zero at the trailing edge. 
In reference 12 this is also shown for Reynolds number equal to 10 a . 
When this assumption was inserted in "the differential, equation for 
$ it was possible to integrate the equation" directly and the 
results obtained by the direct integration were in good agreement 
with, values obtained by the step— by— step integration of the original 
equation. It is believed that for practical applications the 
simplified procedure based on the assumption of a linear variation 
of cf[2it(r/l}l will be sufficiently precise. However, if a body 
of revolution, ■ with shape markedly different over the rear portion 
from that of the Akron, Is to be Investigated for a r&nge of 
transition points and Reynolds numbers it is suggested that repre- 
sentative cases be computed by using a step— by— step Integration 
over the last 20 percent of the chord. Only in this manner can 
the exactness of the assumption concerning the linear variation 
of Cf[2rt(r/l}] be tested. 
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General Remarks 

In certain cases It Is particularly desirable to know the 
variation of either the thickness or the displacement thickness of 
the turbulent boundary layer over a curved surface or a body of 
revolution. In order to determine this precisely from a knowledge 
of the amentum thickness dl strlbut. 1 on It would be necessary., of 
course* to have added information about the variation of the shape 
factor E and the velocity distribution through the boundary layer « 
In reference 24* von Boenhoff and Ttetervin have discussed certain 
aspects of this problem thoroughly* and have presented an empirical 
differential equation that* vheri used ^ vitff'tEo^moTOntuni equation 
and the skin— friction relation*- penmits tracing the development of 
the turbulent boundary layer to the separation point. The calcu- 
lations necessary for the solution of these equations are* however* 
of considerable length. 

The relations in reference 24 make it possible to calculate 
the variation of momentum thickness and boimdary— layer-'Shape factor 
accurately. For airfoils at low speeds* this calculation Involves 
the use of three equations: 



(H+2) 


e au _ 1 
u " pr 


£= 2.555 In (Ik 075 B G f Ff ' 

© m =, © 4 * e ®° ( H - 2.975} T ^ 0 djJ ^ _e. 035 (B-l. 286)1 
dx ■■ ' L U dx J 


The first two of these equations are the basic equations of the 
present report and the last equation is an empirical relation which 
was developed by von Boenhoff and Teterv:.n. It so happens that the 
first equation is quite Insensitive to ..the value of H. Thus for 
nonseparated flow tjbe variation of 6 * as computed by means of the 
method of the present report .on the basis "IbT’tBe’S a 

constant value of E, is quite accurate. Once the chordwise distri- 
bution of 6 is found the solution of the empirical equation of 
reference 24. is simplified considerably. If the numerical integra- 
tion of the third equation above gives a value of H in excess of 
1.8 at any chordwise station in the turbulent regime* the imminence 
of turbulent separation can bo suspected and the method of the 
present report cannot be applied behind thatstation. 

The present report treats the turbulent boundary layer in terms 
of the boundary— layer momentum ’thiclaiessV’ SIhcd*““for an unseparated 
turbulent boundary layer a value of H* the ratio of displacement to 
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momentum thickness, of from 1.4 to 1.5 is indicated at low speeds 
and, since the effect of Mach number an E has been assumed to he 
negligible, the value of the displacement thickness follows imme- 
diately. In order to obtain the boundary— layer thickness, consist- 
ency demands the use of the type of velocity distribution which was 
used in the development of the logarithmic relation between the 
skin friction factor £ and the boundary— layer momentum thickness. 
fEowever, this approach involves theoretical difficulties which can 
be circumvented by the following means. It has been previously 
noted that the logarithmic relationship between the local skin- 
friction coefficient and the boundary— layer momentum thickness 
loads to a variation- of local 'skin-frict : on coefficient with 
Reynolds number which is numerically eq.ual to that predicted by the 
power law developed by Falknar in reference 7 . This power law 
relationship between local skin friction coefficient and boundary- 
layer momentum thickness leads to the concilia ion that the variation 
of velocity through the boundary layer is related to the distance 
normal to the surface by the expression 


u = t'Z ' \ 1/5 

U W 


Using this approximation for the velocity profile It Immediately 
follows the houndary-layer thickness is approximately 

8.4(1 + Mil) times the 'boundary-* layer momentum thickness. 


In reference 4 and 20 it is shown, for a wide variety of air- 
foil sections, that the drag coefficients and ‘boundary— laye r— 
thickness distributions calculated from the Squire and Young equations 
agree well with experimental values. - In the results that have been 
given in the present report. It has been shown that for various 
airfoil shapes and for a representative ran^e of Reynolds numbers 
the calculated value of drag coefficient for each case is very close 
to that obtained by previously established methods, including those 
of Squire and Young. Since the method of the present report Is of 
the same accuracy as previous methods and Is more general and easily 
applied, its use In the calculation of drag coefficients and boundary- 
layer distributions is therefore recommended. 
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APPENDIX A 


Symbols 


General Te rms 

1 ~ 

Cf local skinr-frlction coefficient (2T/p 0 U 0 a ) 

D drag of body (per unit length for airfoil section) 
k thermal conductivity 
M Mach number of free stream 
Pr Pranitl ’ s number / \ 

v * J 


dynamic .pres sure of free stream OgPo^o 2 ) 
u local, velocity irieide boundary layer or in wake 
U Q velocity of undisturbed stream 

U local velocity outside boundary layer or at edge of wake 

U nondimensional velocity ratio (U/U Q ) 

7 ratio Of specific heats [('Cp / 'c T ) = i.k] . ■ ■■ 

8 boundary-layer, thickness . - 

£ Bkin~friction factor .( p u U e /T) a 

H 0 . coefficient of, viscosity in free stream , 

IXy. coefficient of viscosity at wall 

p Q density in free stream 

Py density Just outside boundary layer, 

p density inside boundary layer or vake . 
p nondimensional density ratio (pg/p 0 ) ‘ 

t skin friction per unit area ■ 


Airfoil Sect ions 
c ai rf oil chord 

c^ . section profile— drag coefficient (D/-ip 0 U 0 2 c) 
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H 

*c 

w 

X 

X 

7 

7 

z 

6 * 


'boundary— layer shape parameter (&*/d) 
Reynolds number based on chord length f 


£0^ 

Ho 


) 


variably introduced in equation (19) 

distance along airfoil chord 

variable x in nondimens 3 onal terms (r/c) 

distance measured perpendicularly -o airfoil surface or to 
center line of wake 

variable y in nondimens i onal terms (y/c) __ 

variable introduced in equation (12) 

P ^ ""j 

displacement thickness | f Q ^1 — dy f 


6 

momentum thickness! 

PU_ / 

V • * . 

1 — Ji^dy 


i/° 

PyCr \ 

■*) J 


9 momentum, thickness in nondimens i onal terms (6/c) 
Bodies of revolution 

A surface area of body - - 

total drag coefficient (B/-ip 0 tJ 0 ^A) 

Cfc total drag coefficient (D/-|p 0 U 0 e 7 2/S ) 


H 

l 

r 

r 

V 

v 

x 


' g^O O 

boundary— layer parameter 
length of body 

radius of cross section of body 
radius in nondime ns i onal terms (r l) 

— — — —J 

volume of body 

variable introduced in equation (*1-2) 

distance measured parallel to axiB of body from- stagnation 
point / . 


x variable x in nondimens! onal terms £x/z) 


I 
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distance measured perpendicular to surface or to center line 
of valce 

jr distance y in nondimensionaJ. terms (y /l) 
z variable introduced in equation (37) 
a angle between tangent to generator and axiB of body 

. displacement area of boundary layer 

*• { [ 4 s 0 - ^)( r + 7 008 a H 

^ di splacemsnt area of wake^j ^1 — ^~ ) l r 

B variable of integration .( $ ; 2irr) 

0 nondimenslonal variable of integration {Q ll) ^ . 


$ 


f momentum area of boundary layer 


[" 2n f S ~£H_ (l - £~)(r + y cos a^dy] 
L ° P 0 u o ^ U 0 / V / J 


L ° p<AA 

u 0 /v- 


s J 

i momentum area of wake 

2ir f 6 

— SL. ( 1 

tT’) 7 djr 


L o 

Pq^O ^ 

U 0 / J 


S> 


momentum, area in nondimens i onal terms 



Subscripts -• - .. _ . T _ 

S conditions at . stagnation point 

T.E. conditions at trailing edge of airfoil or tail of body of 
revolution 

T.P. conditions at transition point 
v condition at wall or surface 

o free—stream conditions - - ■ - 
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± conditions in laminar "boundary layer at arbitrary point Xi 

0.8 conditions at SO— percent point on axis of body of revolution 
2 conditions in wake where static pressure is that of free stream 
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■ ■ 7 ." - APESNpiX B 

Computation Procedure for Drag Calculation 
1 . Airfoil Section at Arbitrary Lift Coefficient 


A. Prom known airfoil thickness d * atribution or pressure 

' distribution determine velocity, distribution cor- 
responding to desired lift coefficient. 

B. Estimate transition-point location on each surface by 

' the following: 

(a) For maximum velocity in vicinity of 

leading edge, transition occurs 
at chord wise station correspond- 
ing to velocity decrease of from 
5 to 10 percent of maximum velocity. 

(b) For maximum velocity in vicinity of mid- 

chord’. transition occurs near maximum 
’ ^ “ velocity point or, for large Eeynolds 

: numbers, ahead of maximum velocity 
. ...... . pofnt at chor&wiee station where local 

- . ' boundary— layer Reynolds number attains 

a value of about 8000. 


V - i 1 - 0.35 0 T . p: 


TT ■ T.1T i I 

U T.P. *1 . . L 


f' £ T.P. 

, 4 U s * l7 dx 


. C. From stagnation point to transition point the flow in 
boundary layer is laminar] thus at transition point 

; ! P* T.P. 

(6p)m p 2 a li + 0.26 M 2 (l -O.92 tin p S PJ if ' 17 <& 

l ’ e ’ TTm — s.iv ; i u o 


'T.P. 


)| Qrve __ 

D. Transform variable to wm p. = — - — — — ^-(Qp)™ p ITm p 

1+0 . 152M 2 i * r * 

E. From figure 2 obtain z^.P. == I.60U wj^p^n^ip p^ 

F. For turbulent region, from transition point to trailing 

edge, the growth of z is given by 


*T.E. = Z T.P • 


— P T ’ E '5 , u 4 dx 

[l+0,152Ma]U T . E> 3 4j, >p> 
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where 


P «* 


i 

+ 2=i M 2 (l-tf 2 ) J 7 “ « (1-U 2 )3 


G. At trailing edge obtain Wg^jg corresponding to 

Z T.E. figur® H ? . ? 



Section drag coefficient is obtained from 


c a - 


2(1 + 0.15^ 2 )w t . e> /s \2«a 

— p -(%.s.j 


2. Bodies of Revolution at Zero Angle of Attack 



A. From, body thickness di stri b'u+ ' . in ■ or pressure distri- 

bution determine velocity distribution. 

B. Transition from laminar to turbulent flow occurs 

roughly at chordwiae station where local boundary— 
• • layer EeynoldB number attains a value of 8000. 


U. 


5-3 B Z 

7.17 ; — 


T.P 


r T.Pr 


1-0 . 35M 2 ) 1-1 , 67% . p? . 


ll p T T.P. 

\j fSg®- 


17 d3c 


-i o 


From stagnation point to trans.lt ion point the flow* in 
"boundary ‘ layer is laminar | tiius" at transition point 


(9 p) 


T.P. 


2 = 


da. 


■p TT 9«l^rr 2 

r i u t.p. ^.p. 


I |l + 0.26 M 2 f 1-0.92 UT.P. 2 )} 


r a U s * 17 d^ ' 


D. Transform variable to 


« . (5o) t p .% ; 

.1+0.152M 2 


E.’ From figure 2 obtain zgrp^ *= 1.60*)- Wrp p^In^^p^ 


F. For turbulent region from transition point to 80 - 

percent axial station, the growth of z is given 
by 
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where 


f u T.pA a ^ r T.P.V ,al ^ 

z o.e = ^r.pA *r J V= ) + 

tJo.e r o.s 


1.21 E z 


(l+0.152M 2 )U o . s s 


^T.P, 0,8 


2l 

dx 


■l 

7-1 


p = Cl + •— M s ( 1 -TJ 2 )] »[!+ \ M a U--^ 2 )] 


G-. From figure 2 obtain corresponding to z o>Q 

H. Transform variable to 


4>, 


o.a 


2n r o .8 e o,8 


greo-<,^.s(3-+0-15g 
4.075 SjP 0 .s D o.s 


I. From SO^percent axial station to trailing edge, growth 
of ® is found from 



/‘Uo .8 N 3 * 5 / 

^Pb.e V 

.E/ ' 

" PT.E . * 


0 . 766 ( 2 flr o<8 ) Cb.e'tt o„s a 
ln a (w 0>Q ) PT.E. u T.I. a 


/ T 'Vs) u~y- 

'i.a W.l :/ 


5 dx 


J. The drag coefficient of the body of revolution, based on 
the volume to the two-thirds power, is then 


°D = 


g PT.E. < * > T,-E. 

v a/3 /lS 



3.2 
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TABLE I.- COMPARISON OP SECTION DRAG COEFFICIENT FOR SEPARATE 
AIRFOIL SURFACES COMPUTED BY VARIOUS METHODS 


Reynolds . ! 
number I 



Upper Lower 

Burfaoe | surface 

x/c of T, P, | x/c of T, P. 


1C 6 

10® 

10 ® 

107 

10 7 

10 7 

5 x io 7 


.177 
.376 . 
.177 
.376 

.196 

.396 

*196. . 
.396 


...Section ^gjjoefficieat^cd 

Reference 2 i Reference 4 
Upper j Lower t Upper j Lower 
BUrface j surface surface | surface 


0.0046 

.0041 

.0036 

.0030 

.0026 

.0021 

.0024 

.OO65 

.0052 

,oo4i 

.0031 

.0091 

.0067 

.0057- 

.0038 


0.0046 

.0041 

i .0036 : 
; .0030 ; 

i .0026 [ 
.0021 ! 
.0024 
.0050 
.0041 

.0031 j 

,0023 j 
.0066 j 
.0050 ! 
.,0041 i 

.00 29 


o.oo46 

.0042 

.0034 
.0030 
.0025 
.0021 
.0023 
. .0063 
.0050 
.0040 
.0030 
.0091 
.0065 : 
. 0060 . 
.0038: 


o.oo46 
.0042 
.0034 
.0030 
I .0025 
I .0021 
i .0023 
; .0050 

I .0040 

i .0031 
i .0023 
.0065 
i .0049 
j .oo4i 
.0029 


Present report 

Upper ! Lower 
aurfa.ce ! erurfac© 

1 

1 

1 

0.0047 i 0.0047 
.oo4o ! .oo4o 
.0035 ; .0035 
.0030 j .0030 
.0026 .0026 

.0021 ',0021 

.0024 .0024 

.0063 .0049 

.0051 j .0040 
,oo4i 1 .0031 
.0031 • .0023 
.0089 [ .0064 
.0066 j .0049 
.0059 j - 001 * 2 
.0039 i .0029 


e* 


9osiv *°jt m ravK 
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TABLE II.- COMPARISON OF EXPERIMENTAL DRAG COEFFICIENTS FOE 
NACA 0012 AIEFOIL SECTION AT ZERO ANGLE OF ATTACK 
WITH VALUES CALCULATED BY SEVERAL METHODS 



Section drag coefficient. 

cd __ . 


Reynolds 
number 
X 10“ 6 

Experimental 

measurement 

Calculated 
reference 2 

Calculated 
reference 3 

Calculated, 
reference 5 

Calculated 

present 

report 

2.675 

- 1 

O.OO71 

0.0074 

O.OO67 

0.0069 

i - 

0.0067 

3.780 

.0070 

.0072'" 

•:oo 69 

1 

.fore * 

5.350 

.0068 

.0071 

.00.69 

.0070 

.0068 

7.560 

.0067 

.0071 

.0069 

.0069 

.0067 


■No- 


table IH.- THEORETICAL DRAG COEFFICIENTS OF AKRON 
AIRSHIP SHAPE CALCULATED AT SEVERAL REYNOLDS 
. NUMBERS AND AT VARIOUS TRANSITION— 

• - POINT LOCATIONS 


1 

\ 

< 

Reynolds ■] 
number j 

. . : ..... 1 

h* 

Transition ! 
percent Z 1 

Drag coefficient, Ca i 

Present 

report 

i — 

1 Reference 12 | 

10® 1 

4.6' " 

i 

0.00502 

j 0.00508 1 

io t | 

4.6 

! 

,00343 

• .00335 

10 8 t 

4.6 

i 

i 

.00245 

•! .00235 

10® ' i 

25.7 

i 

.00438 

1 . 00446 

107 •: 

25-7 

t 

i 

.00284 

j - .00279 

10 8 ! 

25-7 

\ 

i 

.00198 

! .00189 

10® ; 

53. ^ 


.00312 

! .00316 

10 7 1 

53.4 


.00173 

• .00176 

10 s j 

53.4 

i 

i 

.00114 

; .00115 


r 

1 
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Figs. 1,3 
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Figure 1.- Airfoil section, with boundary layer. 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTIC8 



Figure 3.- Body of revolution with, boundary layer. 
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Fig. 2« 
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(a) Values of z from. 10® to 10® 


Figure 2 (a to o) Qrapb for evaluating w from z for turbulent boun- 
dary layer* . 
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Fig. 2b 
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(b) Values of z from 10® to 10 7 
Figure 2.- (Continued). 
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Figure 3.- 


(o) Values of z from 10 7 to 10^ 
(Concluded) . 



